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We review our recent work on Chiral Lagrangians out of thermal equilibrium, 
which are introduced to analyse the pion gas formed after a Relativistic Heavy Ion 
Collision. Chiral Perturbation Theory is extended by letting f n be time dependent 
and allows to describe explosive production of pions in parametric resonance. This 
mechanism could be relevant if hadronization occurs at the chiral phase transition. 



1 Introduction 

In the last few years there has been a growing interest in nonequilibrium Quan- 
tum Field Theory, motivated by the experiments on Relativistic Heavy Ion Col- 
lisions seeking the Quark-Gluon Plasma (QGP) and its properties. After the 
promising results obtained by SPS at CERN the torch has been passed on to 
RHIC at the BNL. Roughly speaking, these experiments observe the final par- 
ticle spectra, trying to extract information about the possible QGP formation 
and its subsequent expansion until freeze-out, during which the chiral phase 
transition and hadronization take place. Apart from equilibrium properties 
such as in-medium masses and decay widths or the phase diagram, nonequi- 
librium aspects of the expansion may also be relevant. An important example 
is the hadronization process, which is not yet fully understood. One of the 
scenarios suggested to explain the production of hadrons out of the expanding 
plasma is that of Disoriented Chiral Condensates, where pions develop strong 
fluctuations after the chiral phase transition Era. In this context, two mecha- 
nisms can give rise to pion production. The first is spinodal decomposition, 
which takes place in an early stage after the transition, when it is reasonable 
to assume an initial supercooled state. Thus-Jnng wavelength modes grow 
exponentially yielding strong pion correlations EHj'u. The second is parametric 
resonance, where the oscillations of the a mode in a later stage yield explosive 
pion production. The resulting reheating process yields a. final temperature 
Tf ~ 135 MeV compatible with the freeze-out of hadronsla. Both approaches 
are complementary and have been .studied numerically in great detail in the 
large- A'' limit of the 0(N) model Here we will discuss the extension of 
Chiral Perturbation Theory (ChPT) to include noncquiliLrium effects and its 
application to pion production in parametric resonancel3'B. 
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2 Nonequilibrium effective chiral models 

An effective model to describe QCD well below the chiral scale A x ~ 1.2 GeV 
must be based on the chiral Spontaneous Symmetry Breaking (SSB) pattern 
SU L (N f ) x SU R {N f ) — > SU v (N f ) where N f is the number of light flavours. 
One example is the O(N) model, which reproduces the chiral SSB only for 
N — 4 (Nf — 2) and, due to the strong interaction, is nonperturbative in the 
coupling constant, so that nonperturbative schemes such|-asJarge N must be 
implemented. An alternative for low-energy QCD is ChPTtllilJ which provides a 
well-defined expansion in p/A x where p is any meson energy. Such expansion is 
renormalizable order by order. For instance, the 0(p 2 ) lagrangian, the familiar 
Non-Linear Sigma Model (NLSM), generates 0(p A ) loops whose divergences 
are absorbed by adding to the NLSM a 0(p A ) lagrangian including all possible 
terms compatible with the symmetries, with undetermined coefficients, the 
so-called low-energy constants (LEC), which can be fixed with experimental 
information. Regarding nonequilibrium, the ChPT scheme is best suited for a 
stage where the system is into the broken phase. In that case, it provides us 
with the advantages commented above. For instance, it can be used both for 
Nf = 2, 3. Therefore, as far as pion production is concerned, it will be useful 
in the parametric resonance regime, whereas the O(N) model is more adequate 
to describe the spinodal phase. The first step is then to extend the NLSM out 
of equilibrium. The approach we will follow is to perturb the system from an 
equilibrium state and analyze its subsequent evolution. Since the form of the 
lagrangian is dictated by the chiral symmetry, assuming spatial homogeneity 
and isotropy, the only way to introduce the perturbation.-jto leading order is 
by letting the pion decay constant f v be time dependent 0. Thus, the 0(p 2 ) 
NLSM nonequilibrium action becomes 



where J" indicates that the time integration runs along a Schwinger-Keldysh 
contour in the complex plane, which is the standard technique to formulate 
the non-equilibrium path integral. We fix the initial time to t — 0, where the 
system is in thermal equilibrium at temperature TJj. Thus, f(t<0) — f where 
/ = U(l + 0(p 2 )) and U ^ 93 MeV. We restrict to N f = 2 so that the SU(2) 
matrix U= [(f^-Tr^^I+if-n}/ f(t) where ir a (x, t) are the three pion fields. 
Note also that we work in the chiral limit, i.e, the u, d masses vanish exactly, 
so that we have not included explicit symmetry breaking mass terms. One can 
then proceed with ChPT as usual, expanding the action (Q) in pion fields, which 
generates all possible interaction vertices. The chiral power counting remains 




(1) 
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consistent if /// ~ 0(p), f / f ~ 0(p 2 ) and so on. In this sense, ChPT is more 
appropriate to describe a situation not far from equilibrium such as parametric 
resonance. As explained before, from the moment we consider loop diagrams 
the ChPT scheme requires the 0(p 4 ) lagrangian, whose nonequilibrium version 
is more involved to construct than the 0(p 2 ) NLSM in ([!]). For that purpose, 
it is useful to view the action (Q) as an equilibrium NLSM in a spatially flat 
Robertson- Walker (RW) metric background whose scale factor is just f(t)/f. 
This is easily achieved by rescaling the pion fields to 7r a ///(i) and allows to 
write the lagrangian to any order by rising and lowering indices with the RW 
metric and including all new tensor structures, preserving the symmetry. The 
0(p 4 ) lagrangian built in this way is given intltl In addition to the equilibrium 
lagrangian, it involves two new teoas and hence two new LEC, which have 
been determined experimentally in til. These two terms are crucial since they 
allow to renormalize new time-dependent infinities that were not present in the 
equilibrium case. For instance, the ChPT one-loop diagrams include tadpoles 
proportional to the equal time correlator J ^^d k i Go(t, t, k), where d is the 
space-time dimension and Go(t,t',k) — {Tc^k{t)^k{t')) is the pion two-point 
function. This quantity has a time-dependent divergent part which vanishes 
in equilibrium. Using dimensional regularization, these divergences can be 
absorbed in the new 0(p ) parameters so that observables become finite and 
scale independent 



3 Parametric Resonance in ChPT 

In the nonequilibrium 0(4) model, (a) is time dependent and in the last stage of 
the nonequilibrium evolution it oscillates around the constant vacuum. Thus, 
a(x,t) = &o(t) +5cr(x,i) where 00 (i) is a time-dependent classical background 
which oscillates around ctq = / and 6a includes quantum corrections, which 
are subleading. The oscillations transfer energy to the pion fields, yielding 
exponential growth of the pion correlator via parametric resonance u. The 
connection with ChPT is established through the observation that the NLSM 
corresponds, to leading order, to the 0(4) model with the constraint a 2 + ir 2 = 
a 2 (t) and f(t) ~ ao{t). Therefore, we will take 

/(£) = / [l- | sin Aft] (t>0) (2) 

It is important to bear in mind that using (^) is valid for times t < tsR 
where we denote by tsR the back-reaction time such as the pion correlations 
become typically (7r 2 )(t_Bfl) ~ f 2 ■ When that time scale is reached, (||) is 
no longer a solution of the equations of motion in the 0(4) model and loop 
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Figure 1: The left curve is A uns t(t} / where the dashed line is the tree level contribution 
in (0). The right curve shows /^(t)//7r(0) for the same parameters. 



corrections in ChPT become of the same order as the tree level. It must he 
pointed out that for t ~ tsR, dissipation makes particle production stopp. 
Therefore, the ChPT approach is able to account for nearly all the pion pro- 
duction before tsn- Another point with is worth stressing is that to one loop 
in ChPT, the nonequilibrium chiral power counting demands qM 2 — 0(p 2 ). 
That is, we are studying small oscillations, which in inflationary Cosmology is 
called narrow resonance. The mass parameter M does not need to be small 
and in fact numerically M ~ m ai although there is no need to invoke the a 
particle in the ChPT context. Once the assumptions and limitations of our 
approach are clear, let us sketch how the parametric resonance works in ChPT. 
To 0(ir 2 ), the action (|l|) yields the "free" lagrangian, which now contains a 
time-dependent mass term, so that, taking into account (Q), Go satisfies 

{d 2 + k 2 - (qM 2 /2) sin Aft} G (t, t' , k) = -S c (t - t') (3) 

The above differential equation has the form of the Mathieu equation, 
whose solutions develops instability bands in momentum space, centered at 
k n ~ nM/2 and of width Afc„ = 0(q n ). Therefore, to the order we are 
considering, i.e, 0(q), only the first band is resonant. 

The first observable one can calculate in one-loop ChPT is whose 
NLO corrections include tree level diagrams coming from the 0(p A ) lagrangian 
and tadpoles containing the exponentially growing function 

A unst {t) = ^[ dkk 2 [G (t,t,k)-G o (0,0,k)] (4) 

We have plotted A unst (t) and f%(t) in Figure Q, where-/* ,t (t) correspond 
to the spatial and time components of the axial current 13113. The parameters 
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Figure 2: The left solid curve is the time average of n(k, t) from t = to tsR and the dashed 
line is that of (k 2 /M 2 )n(k,t). The right curve is (n(t))/V = (1/2tt 2 ) J dkk 2 n(k,t). The 
parameters M, q, Ti are the same as in Figure hi 



have been chosen so that tsn. — 10 fm/c, the typical plasma lifetime. We 
have estimated tsn as the time when A„„ si equals the tree level contribution. 
Observe that the f^(t) oscillations are not damped, which is a reminder that 
we have not included the back-reaction properly and should not worry us as far 
as pion production is concerned. We also see that the central value decreases 
with time. This effect can be interpreted as a reheating of the system (since /„• 
decreases with, the temperature) which here leads to a final temperature Tf ~ 
125-140 MeVEl, compatible with hadronic freeze-out. 

The pion distribution function n(k,t) and hence the pion number may 
also be analysed in this context. First, one needs to provide an appropri- 
ate definition for n(k, t) at nonequilibrium since, similarly to QFT in curved 
space-time, the vacuum state is time-dependent. A consistent approach is to 
define it in terms of the energy-momentum tensor expectation value, using a 
suitable point-splitting technique. In this way, one can calculate the number 
of massless pions to one loop in ChPTtl After a careful evaluation of all the 
NLO contributions it turns out n NLO (k,t) = n LO (k,t) + 0(q 2 ). In Figure | 
we have plotted the time average of n(k, t) as well as the total pion density as 
a function of time. The pion spectral function is peaked around the resonant 
frequency k ~ 300 MeV, which could indeed be of importance for RHICEl 
Note that the pion density grows without stop because the back-reaction is 
neglected. However, as commented before, for practical purposes we only need 
the number of pions produced by the end of the expansinn__for which we get 
0.2 per fm 3 . Our predictions agree with the 0(4) modelBu'Ef . 
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4 Conclusions and Outlook 



Chiral Effective Models can be used to describe noncquilibrium phenomena 
such as particle production, which are relevant for the hadronization of a QGP. 
We have showed our results using ChPT, which provides a systematic perturba- 
tive expansion in energies and time derivatives. Explosive production of pions 
takes place when f^if) oscillates in the parametric resonance regime. The 
final pion spectra and pion density agree with recent determinations within 
the 0(4) model. We believe that our results are promising and provide useful 
techniques for future theoretical analysis related to RHIC applications. 
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